Plasma surface dynamics and smoothing in the relativistic few-cycle regime 
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In laser-plasma interactions it is widely accepted that a non-uniform interaction surface will 
invariably seed hydrodynamic instabilities and a growth in the amplitude of the initial modulation. 
Recent experimental results [Dromey, Nat. Phys. 2009] have demonstrated that there must be target 
smoothing in femtosecond timescale relativistic interactions, contrary to prevailing expectation. In 
this paper we develop a theoretical description of the physical process that underlies this novel 
phenomena. We show that the surface dynamics in the few-cycle relativistic regime is dominated by 
the coherent electron motion resulting in a smoothing of the electron surface. This stabilization of 
plasma surfaces is unique in laser-plasma interactions and demonstrates that dynamics in the few- 
cycle regime differ fundamentally from the longer pulse regimes. This has important consequences 
for applications such as radiation pressure acceleration of protons and ions and harmonic generation 
from relativist ically oscillating surfaces. 
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The dynamics of the plasma- vacuum surface are of crit- 
ical importance to the understanding laser-plasma inter- 
actions. Until now, the prevailing expectation has been 
that laser-plasma surfaces are inherently unstable enti- 
ties, since the morphology of the laser-plasma interaction 
surface is determined by hydrodynamics for timescales 
where ion motion is significant and initial perturbation 
will grow rapidly due to instabilities. Specifically, the sur- 
face dynamics play a central role for many applications of 
intense laser-plasma interactions. The best known exam- 
ple is probably the Ray leigh- Taylor instability which lim- 
its the parameter space that can practically be accessed 
in Inertial Confinement Fusion [l| . For relativistic pulses 
with durations of >100fs hydrodynamic motion is also 
the key ingredient to determining the surface morphol- 
ogy. In this regime the role of the light-fluid is performed 
by the ponderomotive pressure of an intense laser @, Q , 
in close analogy to the Ray leigh- Taylor instability. The 
growth and evolution of this instability is essential to un- 
derstanding hole-boring 0, [H| . The sole exception to this 
view of the interaction-surface dynamics were ultrashort, 
few-cycle interactions at modest intensities where there 
is insufficient time for ion motion to take place and con- 
sequently the initial surface shape is typically assumed 
to remain unchanged during the interaction. 

However, recent results by Dromey et al. showed 
that the reflection of high-order harmonics from initially 
rough targets was consistent with the existence of a very 
effective smoothing mechanism 0, 0|. This suggests a 
new paradigm in laser plasma interactions, in that the ex- 
act opposite of the generally accepted behaviour occurs — 
surface smoothing rather than a modulation growth due 
to instabilities. Motivated by this surprising behaviour, 
we have investigated the interaction of a relativistically 
strong laser pulse with an overdense, modulated plasma 



surface. We show for the first time that the surface mor- 
phology in the few-cycle relativistic regime is dominated 
by coherent electron motion in the laser field. This is 
a complete departure from physical picture of how the 
interaction surface evolves in laser-plasma interactions 
and has important consequences for the frontier of ultra- 
fast science in the relativistic regime. This effect exists 
in a broad, practically important parameter range span- 
ning ultrafast pulses with intensities from 10 18 Wcm -2 
upwards. 

To gain insight into the particle motion in the over- 
dense plasma we use a simple one-dimensional model 
[l2l-[l4). This single-particle model describes the motion 
of an incompressible electron layer bound to immobile 
ion background via charge-separation fields under the in- 
fluence of normally incident, linearly polarized electro- 
magnetic wave. This layer, later referred to simply as 
the electron, is initially located on the vacuum-plasma 
interface at x = 0. Taking into account that the charge 
separation fields are proportional to the electron longitu- 
dinal coordinate x, the equations of motion can be readily 
obtained as: 
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where x, y are the propagation and transversal coordi- 
nates respectively, j3 XiV and p XiV - velocity and momenta 
components respectively, a y - the driving vector poten- 
tial and n e is the electron density. We work in relativistic 
units. The normalized quantities for vector potential a, 
time t, length /, momentum p, and density n are obtained 
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from their counterparts in SI- units A, t 1 , l' : p' : and n' via 

eA u L p' n' 
a = , t = ujLt , / = — / , p = , n = . (3) 

m e c c m e c n cr 

Here e and m e are the charge and the mass of the elec- 
tron, ujl is the laser angular frequency, c is the speed of 
light in vacuum, and n cr = s§m e uj\i 'e 2 is the electron 
critical density. 

In equations (j2j) a y (t,x) denotes the driving vector 
potential on the vacuum-plasma interface, which results 
from the interference between the incident and reflected 
wave. It can be found by imposing the standard bound- 
ary conditions for the continuity of electromagnetic fields, 
i.e. of vector potential and its spatial derivative d x a y at 
the plasma- vacuum interface. Without losing the gen- 
erality, the incident a*, reflected a r , and transmitted a 1 
vector potentials can be taken in the form a y (t — x') = 
—Ei • sin(t — x'), a r y (t + x') = —E r • sin(t + x' + (j) r ) 
and a y (t,x') = —E t • sin(t + <p t ) • exp[— uj p {x' — x)] re- 
spectively. With x' we denote the longitudinal coordi- 
nate for the electromagnetic field, while keeping the no- 
tation of x for the coordinate of the electron. Applying 
the boundary conditions (thus setting x' — x) one gets 

Ei = E r = \- yjl + u% • Et and r = 20 £ = 2(a - x), 
where a c± arctanco> p with uj p = ^JrT e the plasma fre- 



quency. One can use the transmitted vector potential at 
x = x' to obtain the driving vector potential: 
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It is important to notice that the actual vector po- 
tential driving the electron is approximately uo p /2 times 
lower than the incident one. As a consequence relativist ic 
effects to the electron motion and corresponding correc- 
tions to the skin depth become important only when the 
amplitude of the incoming light ao exceeds uj p /2 [l2| . 

Results of the model calculations are presented in 
Fig. [TJ The trajectory of an electron interacting with 
a laser pulse having Gaussian envelope of 4-cycles Full 
Width Half Maximum (FWHM) duration tfwhm and 
amplitude ao = 10 (corresponding to an intensity of 
1.37 • 10 20 W/cm 2 for a laser wavelength of = 1 /im) 
is shown in Fig. Plasma density is n e = 400. Fig. [TJd 
shows the transverse coordinate y as a function of time t. 
In Fig.Hh the solid line shows the longitudinal coordinate 
x (horizontal axis) as a function of time t (vertical axis). 

Figure QJ allows us to understand the origin of the 
harmonic generation process. One can see that during 
the interaction the model-electron (the step-like reflect- 
ing surface) oscillates in longitudinal direction with twice 
the laser frequency. Each time the surface moves towards 
the laser it produces a flash with attosecond duration 
[15| . This simple and intuitive picture is called the Oscil- 
lating Mirror model and was proposed by Bulanov et at. 
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FIG. 1: Electron motion obtained using the capacitor model 
for laser pulse with ao = 10 with 4 cycles FWHM-duration 
and n e = 400. Electron is initially located at x e = y e — 0. 
Subflgure (a) shows the electron trajectory, subflgure (b) 
demonstrates the behaviour of the transverse coordinate y e in 
time, on subflgure (c) the dashed line represents the longitu- 
dinal coordinate x e of the electron (vertical axis) versus time 
(horizontal axis) obtained from the model, the color coded im- 
age displays the spatio-temporal picture of the electron den- 
sity obtained from ID-PIC simulations with same laser and 
plasma parameters. 



[8[ and further developed by Lichters et al. [9]. In this 
paper we want to pay attention to the transverse motion 
of the electron, which extends to a considerable fraction 
of the laser wavelength (see Fig. [TJ)) and therefore might 
be responsible for the surface smoothing. Indeed if the 
transverse motion of the electron exceeds the character- 
istic size of the modulations on a rough surface, then the 
roughness is likely to disappear. 

In assessing the role of the transverse motion as a pos- 
sible smoothing mechanism a simple expression for its 
amplitude is needed. Neglecting longitudinal motion one 
can get from Eq. [H an estimate for the amplitude of the 
transverse motion ymax : 
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The dependence of the transverse electron motion am- 
plitude ymax on the laser pulse amplitude ao for plasma 
density n e = 400 is shown on Fig. [2d. The solid line 
shows the results obtained by numerically solving the 
model equations and the dashed line represents equa- 
tion (|5]). The simple estimate (|5]) works fairly good for 
the parameter range studied, and its simplicity makes 
it convenient for the following estimates. More accurate 
results can be obtained by numerically integrating the 
model equations. 

Having estimated the amplitude of the transverse co- 
ordinate y m ax one can establish an ad-hoc criterion for 
surface smoothing to occur based on the ratio of this 
amplitude to the characteristic roughness size h. For in- 
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FIG. 2: Dependance of the amplitudes of longitudinal 
(a) and transverse Umax (b) motion on laser amplitude 
ao. In figure (a) the circles represent the results of ID PIC 
simulations and the solid line shows the results of the numer- 
ical integration of the capacitor model. In figure (b) the solid 
line depicts the numerical integration of the model equations 
and the dashed line is obtained from equation J5}. 



stances where the transverse motion is on the order of the 
characteristic roughness size within the interaction area, 
considerable smoothing can be expected. One can define 
a dimensionless parameter £ separating the case when 
smoothing takes place from the case when the roughness 
survives during the interaction: 
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where h x and h y are the characteristic roughness size in 
longitudinal and transverse directions respectively. We 
make an assumption that the boundary conditions stay 
the same independent of surface structure and that the 
field exponentially decays inside the plasma. In the case 
when £ 1 the roughness according to our criterion 
vanishes. We show further that even in the case when 
£ ~ 1 substantial smoothing is observed. 

In order to check the validity of the afore-described 
model and to demonstrate the surface smoothing we have 
conducted a series of ID and 2D PIC simulations using 
the code PICWIG [16] with clean and rough surfaces for 
different laser amplitudes ao- The code allows the simu- 
lation of the interaction of the intense laser pulses with 
pre-ionized non-collisional plasma with the beam inci- 
dent normally onto the target. The typical plasma den- 
sity used in ID simulations is n e = 400 and n e = 30 
in the 2D case. A step-like vacuum-plasma interface is 
assumed, the ions are immobile. In the 2D case the sur- 
face is modulated sinusoidally in order to simulate the 
roughness (see left part of Fig. |4]). For convenience, the 
modulation period and amplitude are linked and the po- 
sition of the vacuum-plasma interface is given by the law 
x = h • sin(27n//ft). The laser pulse amplitude was varied 
up to ao = 20 in the ID case and is fixed to ao = 10 in 
the 2D scenario. Throughout the paper we use FWHM 
of the electric field as the definition of the laser pulse 
duration and use pulses with an electric field that has a 



Gaussian envelope function in both time and space. 
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where p and tl are the width of the focus and duration 
of the laser pulse respectively. The FWHM duration is 
related to tl by tfwhm = tlV8 ln2. In the ID case the 
size of the simulation box is 7A/,, the time step is T^/1000 
with Tl the period of the driving laser and each plasma 
cell is initially occupied by 1000 macro-electrons. In the 
2D case the size of the simulation box is 3.5Al in laser 
propagation direction and 40 in polarization direction. 
The time step is Tl / 300 and the laser propagation direc- 
tion spatial step is Al/300. Each cell is initially occupied 
by 50 macro-electrons. 

The results of ID simulations are presented on Fig. 
and Fig. [2^. The color-coded image on Fig. presents 
the spatio-temporal evolution of the electron density ob- 
tained from simulations with same laser and plasma pa- 
rameters as in the model (solid line). One can see that the 
model is in agreement with the PIC simulations. Fig. [2^ 
shows results of the model calculations of longitudinal 
electron amplitude (solid line) for different laser ampli- 
tudes ao compared to simulations (circles). The fact 
that simulation results lie on the curve obtained from 
the model and as longitudinal motion is directly corre- 
lated to transverse motion allows us to claim that the 
model works well and gives correct results for both lon- 
gitudinal (Fig. [2^) and transverse (Fig. Et>) coordinates. 
The latter are hard to obtain from ID PIC simulation 
as the particles leave the interaction region and are very 
intricate to trace. 

In the 2D case we investigate the spatial beaming of 
harmonics as a possible indication of smoothing. We an- 
alyze the propagation of the harmonics emission away 
from the target using Kirchhoff diffraction theory 1171 fol- 
lowing the approach used in earlier investigations [Io[ [Hj . 
The harmonic beam (from 15th to 25th harmonic, cen- 
tral wavelength 0.05A/,) 200 Xl away from the target is 
shown on Fig. [3l On all four sub-figures the color surface 
presents the distribution of normalized intensity of the 
filtered harmonics as a function of both time t and trans- 
verse coordinate y (the ceiling panel shows the same data 
as a color-coded image). The upper-right plane shows 
the projection of the beam to the time axis thus the 
time structure of the harmonics beam exhibiting a train 
of several attosecond pulses. On the upper- left plane 
the intensity distribution of the harmonics beam as a 
function of transverse coordinate y is shown (black solid 
line). Results presented on Fig.[3]a,b,c,d are obtained for 
a surface with modulation size h = (smooth surface, 
£ -> oo), h = 0.05A L (£ « 0.6), h = 0.1A L (£ « 0.05) and 
h = 0.2Xl (£ ~ 0.0008) respectively. There are several 
important points to mention. 

First, for the simulation parameters considered the dis- 
tance of 200 Xl corresponds to the position of the har- 
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FIG. 3: Farfield distribution of the reflected harmonics beam 200Al away from the target for a) smooth surface, b) surface 
with modulation size h = 0.05Al, c) surface with modulation size h = O.IAl, d) surface with modulation size h = 0.2Al- 



monies focus due to surface denting as discussed in the 
paper by Horlein et al [10]. This can be illustrated 
from the Fig. [3^i by the fact that the transverse width 
of the reflected harmonics beam (see graph in the upper- 
left plane) is much less than the initial laser width with 
p 5A L . 

Secondly for dimensionless smoothing parameter £ on 
the order of unity the spatial and temporal structure 
of the harmonic beam is not influenced by the surface 
roughness. Figure [3] shows the harmonic orders from the 
15th to the 25th, which should undergo diffuse reflection 
by each of the rough surfaces simulated. Contrary to 
the Rayleigh criterion [ll| , but in agreement with exper- 
imental observation [6| , almost no change in the harmonic 
beam structure is observed for £ w 0.6 (Fig. [3b) in good 
agreement with our ad-hoc smoothing criterion. Surfaces 
with £ <C 1 (Fig. [3fc,d) lead to the speckle-like diffraction 
picture with more energy going to the wings of the beam. 
The fact that the beam is still tolerably collimated hints 
that even though £ <C 1, the characteristic surface rough- 
ness was significantly diminished during the interaction. 
The analysis of the spatial structure of harmonics gener- 
ated on the corrugated surfaces exhibits collimated beam 
structure and serves as an indirect proof of the surface 
smoothing. 
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FIG. 4: (a)Initial density profile and (b) the smoothed density 
profile in the middle of the interaction process for the surface 
with £ 0.6. 

Direct proof of the surface smoothing can be found on 
Fig. E] where initial density distribution (as function of 
longitudinal and transverse coordinates) and the density 
distribution near the moment when the pulse maximum 



reaches the surface are shown (left and right sub-figures 
respectively). The results here are presented for the sur- 
face with £ « 0.6. The evolution of the electron density 
in time can be traced in the animation made from simu- 
lation data (see Supplementary material), showing that 
the transverse motion of the electrons leads to rapid (in 
contrast to the hydrodynamically slow smoothing due to 
ion motion) smoothing of the corrugation. 

In conclusion we have shown for the first time that co- 
herent electron dynamics is the dominant effect shaping 
the laser plasma interaction surface. This is a paradigm 
shift from the way that surface dynamics have been 
viewed to date - as purely hydrodynamic in nature. Due 
to their much smaller mass, electrons can modify the 
surface morphology on the time-scale of even the short- 
est, few-cycle laser pulses and hence must be taken into 
account when considering intense laser-plasma interac- 
tions. This effect has important consequences in the field 
of ultrafast pulses and their application (e.g. harmonic 
generation, ion acceleration via radiation pressure) and 
implies that surface imperfections on a scale smaller than 
the laser wavelength can be neglected. 
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